Radial sensitivity of kaonic atoms and strongly bound K states 
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The strength of the low energy AT^-nucleus real potential has recently received renewed attention 
in view of experimental evidence for the possible existence of strongly bound K~ states. Previous 
fits to kaonic atom data led to either 'shallow' or to 'deep' potentials, where only the former are in 
agreement with chiral approaches but only the latter can produce strongly bound states. Here we 
explore the uncertainties of the J^~-nucleus optical potentials, obtained from fits to kaonic atom 
data, using the functional derivatives of the best-fit values with respect to the potential. We find 
that only the deep type of potential provides information which is applicable to the K~ interaction 
in the nuclear interior. 
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PACS numbers: 13.75.Gx, 21.10.Gv, 25.80.Dj 



I. MOTIVATION AND BACKGROUND 

Since the early days of kaonic atom experiments it was 
known [l|, Q that due to the strength of the nu- 
clear absorption the real part of the iiT^-nucleus poten- 
tial played a secondary role compared to the imaginary 
part of the potential and consequently it could not be 
determined uniquely from fits to the data for a single 
target nucleus. More recent 'global' fits to large sets 
of data encompassing the whole of the periodic table 
showed [3] that although traditional 'ip' potentials yield 
reasonably good fits to the data, the use of phenomeno- 
logical density-dependent t{p) amplitudes leads to signif- 
icantly better fits. When extrapolated into the interior 
of nuclei the real part of the potentials is typically 180 
MeV deep for the density-dependent variety whereas for 
the 'ip' potentials it is typically less than 100 MeV. We 
note that chiral-motivated A'^-nucleus potentials [1, Q 
are shallower than the phenomcnological 'tp^ potentials. 

Figure [1] shows, as an example, the real part of the 
K~ optical potential for Ni obtained from global fits to 
kaonic atoms data using several phenomenological mod- 
els for the interaction. The simplest Hp'' approach where 
the real and imaginary parts of the effective i-matrix are 
determined from fits to the data yields a of 130 for the 
65 data points used in the fit. Adding an adjustable non- 
linear term leads to the deep potential 'DD' of Ref. Q 
with a of 103 and a greatly increased depth. Also 
shown in the figure is another potential using a geomet- 
rical approach to the density-dependence of t (the 'F' 
potential of Ref. f^l), leading to of 84. The similarity 
of the two deep potentials and the great difference com- 
pared to the shallow one are clearly observed. The other 
curve 'FB' with its error band, also of — 84, is dis- 
cussed below. A consequence of the depth of the potential 
is the ability to support a strongly bound state, a ques- 
tion that was highlighted again recently by experimental 
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FIG. 1: Real part of the /S'~-Ni optical potential for various 
models and the values of x^ for 65 data points in parentheses, 
see text. 



reports on candidates for i^-nuclear bound states in the 
range of binding energy ~ 100-200 MeV 0,[i,&[i3l- 
However, very recently a possible explanation of the ex- 
perimental results of FINUDA [To| was given in terms 
of Ap final state interaction [ll| and the observed width 
was shown to be at odds with recent Faddeev calculations 

Obviously, the '<p'-type of potential cannot gener- 
ate strongly bound nuclear states in the energy range 
Bj^ ~ 100 — 200 MeV, whereas the deeper potentials 
('DD' and 'F') might do. In the present work we address 
the question of how well is the real part of the _fi'^-nucleus 
potential determined, with its ability to support strongly 
bound states as the topic of interest. 

Estimating the uncertainties of hadron-nucleus poten- 
tials as function of position is not a simple task. For 
example, in the 'ip' approach the shape of the potential 
is determined by the nuclear density distribution and the 
uncertainty in the strength parameter, as obtained from 
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fits to the data, implies a fixed relative uncertainty 
at all radii, which is, of course, baseless. Details vary 
when more elaborate forms such as 'DD' or 'F' are used, 
but one is left essentially with analytical continuation 
into the nuclear interior of potentials that might be well- 
determined only close to the nuclear surface. 'Model- 
independent' methods have been used in analyses of elas- 
tic scattering data for various particles [13] to alleviate 
this problem. However, applying e.g. the Fourier-Bessel 
(FB) method in global analyses of kaonic atom data end 
up in too few terms in the series, thus making the uncer- 
tainties unrealistic in their dependence on position. This 
is illustrated in Fig. [1] by the 'FB' curve, obtained by 
adding a Fourier-Bessel series to a Hp' potential. Only 
three terms in the series are needed to achieve a of 84 
and the potential becomes deep, in agreement with the 
other two 'deep' solutions. The error band obtained from 
the FB method [l^ is, nevertheless, unrealistic because 
only three FB terms are used. However, an increase in 
the number of terms is found to be unjustified numeri- 
cally in this case. 

In the present work we adopt a functional-derivative 
approach in order to get more realistic position- 
dependence of the uncertainties of ^"-nucleus poten- 
tials. The method is applied to the two types of po- 
tentials ('shallow' and 'deep') mentioned above and it is 
shown that deep potentials within nuclei are reliably ob- 
tained from fits to experimental results for kaonic atoms. 



II. METHOD 

The radial sensitivity of exotic atom data was ad- 
dressed before Q with the help of a 'notch test', intro- 
ducing a local perturbation into the potential and study- 
ing the changes in the fit to the data as function of po- 
sition of the perturbation. The results gave at least a 
semi-quantitative information on what are the radial re- 
gions which are being probed by the various types of 
exotic atoms. In fact, the difference in that respect be- 
tween deep and shallow kaonic atom potentials could be 
observed. However, the extent of the perturbation was 
somewhat arbitrary and in the present work we report 
on extending that approach to a mathematically well- 
defined limit. 

In order to study the radial sensitivity of global fits 
to kaonic atom data, it is instructive to define the ra- 
dial position parameter in a 'natural' way using, e.g., 
the known charge distribution for each nuclear species in 
the data base. The radial position r is then defined as 
r = Rc + rjac, where Rc and are the radius and diffuse- 
ness parameters, respectively, of a two-parameter Fermi 
(2pF) charge distribution . In that way rj becomes the 
relevant radial parameter when handling together data 
for several nuclear species along the periodic table. The 
value of can be regarded now as a functional of a global 
optical potential V{ri), i.e. = X^I^l*?)]- The param- 
eter 77 is a continuous variable, however it is instructive 



to start the discussion of variations by assuming that the 
global optical potential is defined on a discrete set of grid 
points {77^, i — 1 . . . N}, so that depends on a set of 
N parameters Vi — V{rii). The variation of due to a 
small change in these parameters is simply 

The equivalent expression for the continuous function 
V{r]) can be obtained by taking the limit for a very dense 
grid, leading to 

= / '^^^^^^'^^ ' 

where 

Sx'[V{v)] ^ X^[Viv) + eS,{r,-7^')]-x^[V{v)] 

SV{t]') cr^oe^o e 

(3) 

is the functional derivatives (FD) of x^[V^]. The notation 
daifj — rj') stands for an approximated ^-function. From 
Eq. ([2]) it is seen that the FD determines the effect of a 
local change in the optical potential on x^. Conversely 
it can be said that the optical potential sensitivity to the 
experimental data is determined by the magnitude of the 
FD. In practice the calculation of the FD was carried out 
by multiplying the best fit potential by a factor 

! = i + (Mi^-Tj') (4) 

using a normalized Gaussian with a range parameter a 
for the smeared (5-function, 

5.(ry-r;') = -^e-("-'')^/2-^ (5) 
v27rCT 

For finite values of e and cr the FD can be approximated 
by 

5x^[Vm 1 X^[V{ri){l + e5.{rj-r,'))]-x^[V{Tl)] 
5V{rj') ^ V{7j') e 

(6) 

The parameter e was used for a fractional change in the 
potential and the limit e ^ was obtained numerically 
for several values of a and then extrapolated to cr = 
0. Good numerical stability and good convergence were 
obtained in all cases. Here rj and a are dimensionless 
variables. 



III. RESULTS AND CONCLUSIONS 

The _ftr~-nucleus potentials used in the present work 
are taken from recent global fits [6| to kaonic atom data 
from ^Li to U, a total of 65 data points. A two-parameter 
fit with a tp potential yields a total value of 130 
whereas a four-parameter fit with t[p) amplitudes yields 
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FIG. 2: Functional derivatives of x with respect to the full 
complex (C) and real (R) potential as function of ri, where 
r = i?c + rjac, with J?c and the radius and diffuseness 
parameters, respectively, of the charge distribution. Results 
are shown for the tp potential and for the t[p) 'F' potential 
ofRef. Q. 
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FIG. 3: Real potentials (top) and FD (bottom) for the 'F' 
potential (left) and the tp potential (right) for interaction 
with Ni. The regions between the vertical dotted lines indicate 
where the potentials are determined reliably, see text. 



X^=84 for the 65 data points. In calculating FD for 
global fits radial positions were defined in terms of units 
of diffuseness relative to the charge radius, as described 
above. For several of the lightest nuclei harmonic oscilla- 
tor densities are more appropriate and had indeed been 
used in the fits of Ref. Q. These nuclei have been ex- 
cluded from the calculations of FD for the global fits so 
as to have full consistency in the use of 2pF density dis- 
tributions. That left 50 data points in the calculations 
of the FD. 

Figure [2] shows the FD for relative variations in the 
real potential and in the full complex potential. Inspect- 
ing the FD with respect to the imaginary potential (not 
shown) reveals additivity of the FD, hence the differences 
between the FD for the complex potential and for the real 
potential in the figure are the corresponding FD with re- 
spect to the imaginary potential. We note that the FD 
for the deep 'F' potential is dominated by the real part of 
the potential whereas for the 'ip' potential the two parts 
make similar contributions to the FD. The appearance of 
regions with negative FD need not be surprising because 
it means that some local variations in the shape of the po- 
tentials may cause further reduction in the values of x^. 
Such variations are impossible in the tp potential. How- 
ever, such variations are included when the more flexible 
'F' potential is introduced, and the minimisation process 
essentially follows implicitly the various FD. We avoid at 
present making a quantitative use of the local values of 
the FD, rather we identify with the help of the FD the 
radial regions to which the kaonic atom data are sensi- 
tive. 

From Fig. [5] it can be inferred that the sensitive region 
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FIG. 4: Same as fig. |3]but for ^^C. 



for the real tp potential is between rj = —1.5 and -q — & 
whereas for the F potential it is between rj = —3.5 and 
77 = 4. Recall that = —2.2 correspond to 90% of the 
central charge density and = 2.2 correspond to 10% of 
that density. It therefore becomes clear that within the 
tp potential there is no sensitivity to the interior of the 
nucleus whereas with the t{p) 'F' potential, which yields 
greatly improved fit to the data, there is sensitivity to 
regions within the full nuclear density. 
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Figure [3] shows similar results for Ni, this time the 
radial variable is r, the radial position. Here a is in 
units of fm and therefore the FD is in units of fm~^. On 
the left are shown the real potential (top) and the FD 
(bottom) for the real 'F' potential and on the right are 
shown the corresponding quantities for the tp potential. 
The regions between the two vertical dotted lines indicate 
where variation of the potential will affect the fit to the 
data, thus suggesting the regions where the potentials 
are determined by experiment. The vertical dashed line 
indicates 50% of the central charge density. It is evident 
that for the 'F' potential the sensitivity extends to depths 
of 180 MeV at radii where the density is essentially the 
full nuclear density. Very similar results are obtained 
also for heavier nuclei, where, e.g. for Pb the distinction 
between the tp and the 'F' potential is even greater than 
for Ni, with respect to the sensitivity of the experiment 
to depths and densities. 

Finally, in Fig. 2] are shown similar results for ^^C, 
which is one of the targets studied experimentally [Toj 
and theoretically jB]- Again it is seen that with the 'F' po- 
tential the sensitive region is at smaller radii and higher 
densities compared to that for the tp potentials. 

In summary, it is found that the 'deep' X-nucleus po- 
tentials, which yield excellent fits to all the kaonic atom 
data, are determined reliably to depths of 150-180 MeV 



at regions of almost the full nuclear matter density and 
consequently they may be reliably extrapolated to the 
nuclear interior. In contrast the 'shallow' type of po- 
tentials, which yield inferior fits to the data, are well 
determined only at the nuclear surface and one cannot 
infer from these what is the depth of the i^-nucleus po- 
tential in the nuclear interior. The different sensitivi- 
ties result from the potentials themselves: the additional 
attraction provided by the deep potentials enhances the 
atomic wavefunctions within the nucleus Q thus creating 
the sensitivity at smaller radii. We conclude that optical 
potentials derived from the observed strong-interaction 
effects in kaonic atoms are sufficiently deep to support 
strongly-bound antikaon states, but it does not necessar- 
ily imply that such states do exist. Moreover, the dis- 
crepancy between the very shallow chiral-motivated po- 
tentials 0, and the deep phenomenological potentials 
remains an open problem. 
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A simple and convenient introduction to the concept of 
func tional derivative can be found in the wikipedia web 
site: http:/ / en. wikipedia. org/wiki/PunctionaLderivative 



